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Institute of Experimental Physics, Warsaw University, 00-681 Warsaw, Poland 

In this paper we report on the measurement of the rate of inclusive tt" production induced by 
charged- current neutrino interactions in a CsHg target at a mean energy of 1.3 GeV in the K2K near 
detector. Out of a sample of 11,606 charged current neutrino interactions, we select 479 tt" events 
with two reconstructed photons. We find that the cross section for the inclusive tt" production 
relative to the charged-current quasi-elastic cross section is 

-^^^ = 0.426 ± 0.032(sfat.) ± 0.035(sj/sf.) 

C^CCQE 

The energy dependent cross section ratio is also measured. The results are consistent with previous 
experiments for exclusive channels on different targets. 

PACS numbers: 13.85.Ni,13.85.Qk 



I. INTRODUCTION 

After the observation of solar neutrino and atmo- 
spheric neutrino oscillations [H-Q and their confirmation 
respectively at reactors [H and accelerators the pri- 
mary aim of current and future neutrino experiments is 
to measure the ^13 mixing angle, and to improve accuracy 
in the measurement of oscillation parameters. One of the 
largest limitations of accelerator-based neutrino experi- 
ments comes from the poor experimental knowledge of 
neutrino cross sections in the GeV energy range. Con- 
cerning the measurement of (?i3 via sub-leading z^^ — > Vg, 
oscillation searches, one of the main backgrounds to the 
Ve signal comes from z^^ neutral-current (NC) interac- 
tions producing tt'^'s. Experimental input on the rate of 
the related charged-current (CC) channel, which is the 
focus of this paper, and measurement of the tt*^ produc- 
tion momentum spectrum, allows better understanding 
of this background. Concerning the improvements in the 
measurement of oscillation parameters, and in particular 
of the atmospheric mixing angle ^23 and mass-squared 
difference Am^g via the measurement of the distortion 
of the neutrino energy spectrum induced by neutrino os- 
cillations, knowledge of the overall yield and interaction 
type composition of CC inelastic interactions is crucial. 
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This is because the reconstruction of neutrino energy in 
CC interactions via kinematic means is less accurate in 
inelastic interactions, compared to quasi-elastic (CCQE) 
interactions. Charged-current inclusive tt" production 
(CCtt^) constitutes a large component of all CC inelas- 
tic interactions. In addition, since uncertainties in the 
nuclear models play a significant role in the neutrino- 
nucleus cross section, it is important to have measure- 
ments on different target materials. 

Although there are several theoretical approaches to 
model these processes, the experimental constraints are 
rather weak. Very little data exists in the few GeV neu- 
trino energy range. Experimental measurements of neu- 
tral pion production via CC interactions of few-GcV neu- 
trinos on de uterium have been collected in the past for 
single pion 



141 final states. At 



[10|-[13| and two pion 
higher energies, CC single tt" production cross sections 
have been measured on deuterium [l5l | and heavy freon 
[l6| targets. 

In this paper we present the measurement of inclusive 
CC neutrino interactions with a tt*^ in the final state made 
with the K2K SciBar/EC detector system. The measure- 
ment presented here is the first result on a carbon target 
in the few-GeV neutrino energy range, and improves the 
precision of previous results on different targets. First, 
we obtain the cross section for this process with respect 
to the cross section for both CC QE and inelastic interac- 
tions. We quote our result as a cross section ratio rather 
than as an absolute cross section in order to reduce the 
impact of large uncertainties in the estimation of the K2K 
neutrino flux affecting the SciBar detector. Second, by 
reconstructing the neutrino energy of CC interactions re- 
sulting in inclusive tt*^ production, we present the energy 
dependence of this cross section ratio. Third, by using 
previous K2K experimental input on CC single pion pro- 
duction (l7j , we interpret our result as a measurement of 
the CC deep inelastic cross section, relative to the CC QE 
cross section. Fourth, we present relevant tf" production 
kinematic distributions of our CCtt^ candidate events. 

The paper is organized as follows: Section HIl describes 
the experimental setup: the neutrino beam and the neu- 
trino near detector at KEK; Section lllll dcscribes the sim- 
ulation of the experiment, focusing on the neutrino inter- 
action simulation. In Section llVI we discuss the ingrcdi- 
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ents of our main cross section analysis, describing the ex- 
perimental signature, the CC event selection, the photon 
selection, the tt*' mass and the neutrino energy recon- 
struction; Section |V] describes our likelihood fit method; 
Section IVll describes the systematic uncertainties affect- 
ing our measurement; Section IVIII presents the energy- 
dependent and energy independent cross section results, 
and the comparison with the neutrino interaction simula- 
tion and with existing results are given in Section fVlI Al 
Conclusions are given in Section IVlIII 



II. EXPERIMENTAL SETUP 
A. Neutrino Beam 

The KEK to Kamioka (K2K) experiment 0, [TlSl is 
a long baseline neutrino oscillation experiment in which 
a beam of muon neutrinos created at KEK is detected 
250 km away in the Super-Kamiokande detector, located 
in Kamioka, Japan. To produce the neutrino beam, pro- 
tons are accelerated by the KEK proton synchrotron to a 
kinetic energy of 12 GeV and then extracted every 2.2 s 
in a single turn to the neutrino beam line. The duration 
of an extraction, or spill, is 1.1 ^s and each spill con- 
tains 9 bunches of protons at a 125 ns time interval. The 
protons are steered to the neutrino beam line to strike 
an aluminum target, producing secondary particles. Two 
toroidal magnetic horns focus the positively charged par- 
ticles, mainly tt+'s, in the forward direction. The focused 
positive pions are allowed to decay into a 200 m long 
tunnel where they produce neutrinos via n'^ — >■ ^I'^i^^. 
A beam absorber made of iron, concrete, and soil is lo- 
cated at the end of the decay volume to stop all particles 
except neutrinos. The direction and intensity of the neu- 
trino beam are checked spill-by-spill by monitoring the 
muons produced by pion decay. The energy spectrum of 
the neutrino beam is checked by occasionally monitoring 
the pions focused by the horn magnets 

Over the duration of the K2K experiment, a total of 
9.2 X 10^^ protons were delivered to the target to generate 
the neutrino beam. The SciBar and EC detectors took 
data from October 2003 until November 2004; 2.02 xlO^^ 
protons on target were accumulated during this time. 

A Monte Carlo (MC) simulation is used to predict the 
properties of the neutrino beam. According to the sim- 
ulation, the beam at the near detector is about 97.3% 
pure with a mean energy of 1.3 GeV. A fit of data 
of neutrino interactions in all the near detectors is used 
to fine-tune the simulated neutrino energy spectrum 
Figure [T] shows the energy spectrum for all muon neutrino 
interactions in the fiducial volume of the SciBar detector. 



B. Neutrino Detectors at KEK 

The near detector system is located 300 m downstream 
of the proton target. The purpose of the near detector 
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FIG. 1. The energy spectrum for all muon neutrino interac- 
tions in the SciBar fiducial volume. The gray boxes corre- 
spond to the shape systematic uncertainty. 



is to measure the direction, flux and energy spectrum of 
neutrinos at KEK before oscillation. The near detector 
is also used for measurements of neutrino cross sections. 

A schematic view of the near detector is shown in 
Fig. [2j The near detector consists of a one kiloton water 
Cerenkov detector (IKT) [2l|, a scintillating-fiber/ water 
target tracker (SciFi) [22], a fully active scintillator-bar 
tracker (SciBar) complemented by a lead and fibers elec- 
tron catcher (EC) and a muon range detector (MRD). In 
this section we describe the SciBar, EC and MRD since 
data taken from these detectors are used in the present 
analysis. A full description of the K2K near detectors 
can be found in 



SciFi Detector 



SciBar Detector 

EC Detector 
Muon Range Detector 




FIG. 2. Schematic view of the near neutrino detector. 
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1. Scintillator Bar detector (SciBar) 

The SciBar detector acts as a fully active neutrino tar- 
get and its primary role is to reconstruct the neutrino 
interaction vertex and detect the final state charged par- 
ticles. 

SciBar consists of 14,848 extruded scintillator 

bars of 1.3 x 2.5 x 300 cm"^. Groups of 116 bars arc 
arranged horizontally or vertically to make one plane. 
The planes are arranged in 64 layers orthogonal to the 
beam, each consisting of one horizontal and one vertical 
plane. The total volume is 3 m x 3 m x 1.7 m, for a 
total mass of ^15 tons. Figure [3] shows a diagram of the 
SciBar detector. 




FIG. 3. Diagram of SciBar and of the Electron Catcher (EC). 

The extruded scintillator bars are produced by FNAL 
[2^. The bars are made of polystyrene (CsHs), PFO 
(1%), and FOFOF (0.03%). Each bar is 1.3 cm x 2.5 
cm X 300 cm and has a 0.25 mm thick reflective coating 
made of Ti02. The peak of the emission spectrum for the 
scintillator is at 420 nm. A 1.5 mm diameter wavelength 
shifting (WLS) fiber (Kuraray Y11(200)MS) is inserted 
in a 1.8 mm hole in each bar to guide the scintillation 
light to multi-anode photomultiplier tubes (MAFMTs). 
The average attenuation length of the WLS fibers is ap- 
proximately 350 cm. The absorption peak for the fibers 
is at 430 nm (matching the emission peak for the scintil- 
lator), and the emission peak is at 476 nm. The scintil- 
lation light produced is detected by Hamamatsu H8804 
MAPMTs. Each MAPMT has 64 channels arranged in 
an 8x8 array. Each pixel is 2 mm x 2 mm. The cathode 
material is Bialkali, with a quantum efficiency of 21% 
at a wavelength of 390 nm. The cathode is sensitive to 
wavelengths between 300 and 650 nm. A typical chan- 



nel gain is 6x10^ at a supply voltage of 800-900 V. The 
basic properties such as gain and linearity are measured 
for each channel before installation. The non-linearity of 
the output signal vs. input charge is 5% at 200 photo 
electrons (p.e.) at a gain of 5x10^. Crosstalk in the 
MAPMT is approximately 3% in neighboring channels. 
Groups of 64 fibers arc bundled together and glued to an 
attachment to be precisely aligned with the pixels of the 
MAPMT. SciBar's readout system [2^ consists of a front- 
end electronics board (FEB) attached to each MAPMT 
and a back-end VME module. The front-end electronics 
uses VA/TA ASICs. The VA is a 32-channel pre-amplifier 
chip with a shaper and multiplexer. The TA provides 
timing information by taking the "OR" of 32 channels. 
Each FEB uses two VA/TA packages to read 64 analog 
signals and two timing signals for each MAPMT. Each 
back-end VME board controls the readout of eight FEBs. 
Flash ADCs are used to digitize the charge information, 
and TDCs are used to process the timing information. 
The pedestal width is approximately 0.3 p.e., and the 
timing resolution is 1.3 ns. In order to monitor and cor- 
rect for gain drift during operation, SciBar is equipped 
with a gain calibration system using LEDs [13] . The sys- 
tem shows that the gain is stable within 5% for the entire 
period of operation. Cosmic ray data collected between 
beam spills are used to calibrate the light yield of each 
channel. The average light yield per bar is approximately 
20 p.e. for a minimum ionizing particle. The light yield 
is stable within 1% for the whole period of operation. 
Pedestal, LED, and cosmic-ray data are taken simultane- 
ously with beam data. A crosstalk correction is applied 
to both data and MC before event reconstruction [l7j . 
After the crosstalk correction, scintillator strips with a 
pulse height larger than two p.e. (corresponding to about 
0.2 MeV) are selected for tracking. Charged particles are 
reconstructed by looking for track projections in each of 
the two-dimensional views (x-z and y-z) using a cellu- 
lar automaton algorithm |28| . Three-dimensional tracks 
are reconstructed by matching the z-edgcs and timing 
information of the 2D tracks. Reconstructed tracks are 
required to have hits in at least 3 consecutive layers. The 
minimum length of a reconstructible track is, therefore, 8 
cm, which corresponds to a momentum threshold of 450 
MeV/c for protons. The reconstruction efficiency for an 
isolated track longer than 10 cm is 99%. The efficiency 
is lower for multiple track events due to overlapping of 
tracks in one or both views. 



2. Electron Catcher (EC) 

The EC detector is an electromagnetic calorimeter in- 
stalled just downstream of SciBar as shown in Fig. [31 The 
main purpose of the EC is the longitudinal containment 
of the electromagnetic showers since the whole SciBar 
corresponds to only 4 radiation lengths. The EC provides 
11 radiation lengths and has 85% energy containment at 
3 GcV. The EC consists of one plane of 30 horizontal 
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modules and one plane of 32 vertical modules. The two 
planes have a cross sectional area of 2. 7m x 2. 6m and 
2.6 m X 2.5 m, respectively. The modules were originally 
made for the CHORUS neutrino experiment at CERN 
j29| . Each module is a sandwich of lead and scintillating 
fibers, built by piling up extruded sheets of grooved lead 
with scintillating fibers positioned in the grooves. A mod- 
ule consists of a stack of 21 lead sheets, 2650 mm long, 
and 740 fibers of 1 mm diameter and 3050 mm long. The 
groove diameter is 1.1 mm and the sheet thickness is 1.9 
mm. The sheets material is 99% lead with 1% antimony 
content to improve its mechanical properties. The stack 
is kept together by a welded steel case. An overall thick- 
ness non-uniformity of less than 2% was achieved through 
the extrusion process. At both ends fibers are bundled in 
two independent groups, defining two different readout 
cells of about 42 x 42 mm^ transverse cross section. The 
fibers are manufactured by Kuraray (type SCSF81) and 
consist of a polystyrene core surrounded by a 30 /xm thick 
acrylic cladding, with an emission maximum in the blue, 
around 420 nm. To improve the light collection unifor- 
mity, an acrylic black paint is applied on the surface of 
the last 5 cm of fibers on each side. This has the effect of 
reducing the light coming from the cladding, which has a 
smaller attenuation length. In addition, in order to select 
the spectral component with a larger attenuation length, 
a yellow filter (Kodak Wratten No. 3) is used. The atten- 
uation length was measured to be (462 ± 53) cm when 
the modules were built [2^ and was recently measured to 
be (400 ± 12) cm. At both ends of the readout cell, fibers 
are grouped into two bundles of hexagonal cross-section 
(22.2 mm apex to apex) and are coupled to a Plexiglas 
light guide, also with hexagonal cross section (24 mm 
apex to apex). The hexagonal shape and the length of 
the light guide were chosen to reduce disuniformities in 
the mixing of the light coming out of the individual fibers 
[soj . The light guides are coupled to 1-1/8 inch diame- 
ter photomultipliers, type R1355/SM from Hamamatsu, 
with a special green extended photocathode, of 25 mm ef- 
fective diameter. The cathode material is Bialkali, with a 
quantum efficiency of 27% in the wavelength range 350- 
450 nm. The cathode is sensitive to wavelengths from 
300 nm to 650 nm. A typical current amplification is 
2.1 X 10^ at the supply vohage of 1600 V. The anode 
dark current is 10 nA. The PMT gain of each channel 
was measured before installation. The non-linearity of 
the output signal vs input charge is 2% at 60 mA (cor- 
responding to 600 photo electrons) at a gain of 2 x 10^. 
The PMT produces a differential signal using the outputs 
of the cathode and the last dynode. Signals are read via 
multipolar differential screened cables, 100 m long. The 
readout system consists of 8 QDC VME (CAEN V792) 
with 32-channel 12-bit ADC. Impedance matching cards 
(CAEN A992 custom modified) are used to convert the 
differential signals into single ended signals and to decou- 
ple the PMT's and the QDC grounds. Cosmic rays mea- 
sured during normal data taking in between the neutrino 
spills are used to calibrate the detector and to monitor 



the gain stability. After the calibration the spread in the 
individual channel response was stable within 1%. The 
pedestal width is approximately 0.7 photo electrons and 
the energy resolution was measured in a test beam as 
14%/V£^(GeV). The energy deposited is reconstructed 
by searching for clusters of nearby hits above threshold. 
In this analysis clusters are reconstructed searching for 
hits with more than 20 MeV in the vertical plane and 10 
MeV in the horizontal plane. Hits in the nearest coun- 
ters are iteratively added to the cluster if their energy is 
greater than 10 MeV(5 MeV) for the vertical(horizontal) 
plane. The cluster position is the energy weighted av- 
erage of the positions of the counters belonging to the 
cluster. 



3. Muon range detector (MRD) 

The MRD [sH is the most downstream detector. It 
consists of 12 layers of iron between 13 layers of vertical 
and horizontal drift-tubes. Each layer is approximately 
7.6 m X 7.6 m. To have good energy resolution over the 
entire energy spectrum, the four upstream iron layers 
are each 10 cm thick, while the other eight planes are 20 
cm thick. The total iron thickness of 2 m covers muon 
energies up to 2.8 GeV, which corresponds to 95% of all 
the muons produced by neutrino interactions in K2K. 
There are 6,632 aluminum drift tubes filled with PIO gas 
(Ar:CH4 = 90%10%). The total mass of the iron is 864 
tons, and the mass of the drift tubes is 51 tons. The 
MRD is used to monitor the stability of the neutrino 
beam direction, profile and spectrum by measuring the 
energy, angle and production point of muons produced by 
CC neutrino interactions in the iron target. The MRD 
is also used to identify muons produced in the upstream 
detectors. The energy and angle of the muon can be 
measured by the combination of the MRD and the other 
fine-grained detectors. It is necessary to measure the 
muon energy and direction in order to reconstruct the 
energy of the incident neutrino for CC events. The MRD 
tracking efficiency is 66%, 95%, and 97.5% for tracks that 
traverse one, two and three iron layers, respectively; for 
longer tracks, the efficiency approaches 99%. A track 
that hits less than three layers of MRD is called " one layer 
hit" (MRDIL), while a track that hits more than three 
iron layers will be reconstructed as a 2D track in xz or yz 
planes. The 2D-tracks pair which has the longest overlap 
is taken as a 3D-track (MRD3D). The range of a track 
is estimated using the path length of the reconstructed 
track in iron. The muon energy is calculated by the range 
of the track. The uncertainty in the muon energy due to 
differences among various calculations of the relationship 
between muon energy and range is 1.7%. The uncertainty 
in the weight of the iron is 1%. Thus, the systematic error 
in the MRD energy scale is quoted as the sum of these 
uncertainties, 2.7%. The energy resolution is estimated 
by Monte Carlo simulation to be 0.12 GeV for forward- 
going muons. The angular resolution is about 5 degrees. 
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III. SIMULATION 



Neutrino Interactions 



The neutrino interaction simulation plays an impor- 
tant role for estimating the event yields, and the topo- 
logical and kinematical properties, for CC neutrino in- 
teractions in SciBar producing neutral pions, as well as 
for background processes. We use the NEUT program 
library to simulate neutrino interactions with protons 
and carbon nuclei within the SciBar detector material. 
NEUT simulates neutrino interactions over a wide 
energy range, from ■~100 MeV up to TeV neutrino ener- 
gies, and on different nuclear targets. 

In the simulation program, the following CC and 
NC neutrino interactions are considered: QE scattering 
{u N I N'), single pion production {v N ^ I N' to), 
coherent tt production {v ^^C — !■ I it ^^C), and deep in- 
elastic scattering (DIS, vN I N' hadrons). In these 
reactions, N and N' are the nucleons (proton or neutron), 
I is the lepton (either a charged Icpton or a neutrino), 
and m is a meson. If the neutrino interaction occurs in 
a carbon nucleus, the interactions of the generated par- 
ticles with the remaining nucleons of the nucleus are also 
simulated. 

The total charged-current cross section predicted by 
NEUT, together with the QE scattering, single pion pro- 
duction and deep inelastic scattering contributions, are 
shown in Fig. 21 overlaid with data from several experi- 
ments. 

Given the K2K beam neutrino energy spectrum. Ta- 
ble U shows the fraction of interactions in SciBar that are 
expected to be QE, single pion, etc. according to the 
simulation. 



Interaction type 


Percent of Total 


Charged-current (CC) 


71.8% 


Uf_,n — > fj,^p (QE) 


32.2% 




18.0% 


Ui_^n — !> fi^ mr^ 


6.2% 


Uf^n — >■ ^Tpn'^ 


5.0% 


DIS 


6.8% 


others (K and rf) 


3.6% 


Neutral-current (NC) 


28.2% 



TABLE I. Expected neutrino interactions in SciBar. 



1. TT -producing charged-current neutrino interactions 

For the simulation of CC neutrino interactions result- 
ing in inclusive tt'^ production, we adopt distinct models, 
depending on the invariant mass W of the hadronic sys- 
tem in the final state, and on the pion multiplicity. A 
summary of the models used to simulate the cross sec- 
tion and the final state kinematics is given in Table [ill 
while more details are given in the text below. 
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FIG. 4. Charged-current total cross section divided by 
the neutrino energy for neutrino-nucleon charged-current 
interactions [s^j. The solid line shows the calculated total 
cross section. The dashed, dotted and dash-dotted lines show 
the calculated quasi-elastic, single pion and deep inelastic 
scattering, respectively. The data points are taken from the 
following ex£eriment_s^ (A) ANL[33|, (Q) GGM77[3l|, (•) 



GGM79 (a)j35| 
(*) BNL86[iy 
IHEP-JINR96 
and (A) BNL82[4, 



(b)M, (*) Serpukhov[37|, (0) ANL82[Il|l, 
_(■) CCFRgofli, (T) CDHSW87|3il, (x) 
[+) IHEP-ITEP79[43l, (□) CCFRR84[43|, 



w 


Process 


Cross 


Final State 


(GeV/c^) 




Section 


Kinematics 


<2.0 


VfiU — !> fl^p-K^ 


M 


[45], isotropic 


1.3-2.0 




[46, 47] 


[48] 


>2.0 


iy^,N ^.-N'-K^X 


[46, 47] 


[49] 



TABLE II. Models used to simulate the cross section and final 
state kinematics for CC inclusive tt'' production. In the table, 
W stands for the invariant mass of the final state hadronic 
system, and N' for nucleons, tt for at least one charged or 
neutral pion, and X for any meson (including none). See text 
for details. 



For W < 2 GeV, and production of single tt'^'s and no 
other pions (charged or neutral), we use the resonance- 
mediated Rein-Sehgal model [4J]. In this model, the in- 
teraction simulation is performed via a two-step process. 
First, the neutrino- induced excitation of the baryon res- 
onance N* is modeled: 



N* 



which is then followed by the resonance decay to a pion- 
nucleon final state: 



N* 



tt" +p. 
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The same 7r°-p final state can be fed by several reso- 
nances. All baryon resonances with W < 2 GeV/c^ are 
taken into account with their corresponding resonance 
width and including possible interferences among them. 
Single K and ij productions are simulated by using the 
same framework as for the dominant single tt production 
processes. The model contains a phenomcnological pa- 
rameter (the single pion axial vector mass, Ma), that 
must be determined experimentally. As the value of Ma 
increases, interactions with higher values (and there- 
fore larger scattering angles) are enhanced. The Ma pa- 
rameter in our resonance-mediated Rein-Sehgal model is 
set to 1.1 GeV/c^. To determine the final state kinemat- 
ics in the decay of the dominant resonance ^33(1232), 
Rein's method [i^ is used to generate the pion angu- 
lar distribution in the resonance rest frame. For the 
other resonances, the directional distribution of the gen- 
erated pion is set to be isotropic in the resonance rest 
frame. The pion angular distribution for the related 
lyp pn'^ mode has been measured [T^], and the re- 

sults agree well with our model. To describe nucleons 
bound in carbon nuclei, nucleons are treated as quasi- 
free particles in motion using a relativistic Fermi gas 
model [5^ [m with 225 MeV/c Fermi surface momen- 
tum, and assuming a nuclear binding energy of 27 MeV. 
The Pauli blocking effect in the decay of the baryon res- 
onance is taken into account by requiring that the mo- 
mentum of the nucleon should be larger than the Fermi 
surface momentum. In addition, pion-less decay for the 
dominant A resonance (A7V — ;> NN) is considered. In 
this case, which is expected to occur with a 20% proba- 
bility, no pion is present in the final state; only a lepton 
and a nucleon are emitted js^ . 

For the production of tt^'s in association with other 
pions (charged or neutral), the deep inelastic scatter- 
ing cross section formalism combined with GRV94 par- 
ton distribution functions j46| is used. Additionally, we 
have included the corrections in the small region de- 
veloped by Bodek and Yang [l^l- For the simulation 
of DIS final state kinematics in which the hadronic in- 
variant mass, W, is larger than 2 GeV/c^, wc use the 
PYTHIA/JETSET library [liL For W <2 GeV/c^, we 
use a custom made program |48| . based on data-driven 
average pion multiplicities and KNO scaling. In the lat- 
ter case, the multiplicity of pions is required to be larger 
than one because single pion production is already taken 
into account by the resonance-mediated single pion pro- 
duction. 



2. Other neutrino interactions 

Resonance-mediated single pion production and deep 
inelastic scattering CC processes that do not result in 
the production of tt^'s, but possibly in the production of 
other mesons such as tt^ , arc also simulated according to 
the models described in Scc. lIII AD The same models are 
used as well to simulate the corresponding NC channels. 



The formalism of CC and NC QE scattering off free 
nucleons used in the simulation is described by Llewellyn- 
Smith [5^ . There is only one parameter in the model to 
be determined experimentally, the QE axial vector mass. 
Ma- As for single pion production via baryon resonances, 
Ma is set to 1.1 GeV/c^ in our simulation, based on near 
detector data p^ . 

Coherent single pion production, that is the interac- 
tion between a neutrino and the entire carbon nucleus 
resulting in the production of single pions and no nuclear 
break-up. is simulated using the formalism developed by 
Rein and Sehgal 5J]. The coherent pion production ax- 
ial vector mass is set to 1.0 GeV/c^ in our model. Only 
neutral-current coherent pion production interactions are 
considered, because the cross section of the CC coherent 
pion production was found to be very small at K2K beam 
energies [ssj . 



3. Intra-nuclear hadronic interactions 

The intra-nuclear interactions of the mesons and nucle- 
ons produced in neutrino interactions with carbon nuclei 
arc also important for this analysis. Due to the propa- 
gation in the nuclear matter of the target nucleus, the 
final state particles observed differ from the one pro- 
duced at the weak interaction vertex. Particles absorp- 
tion or production as well as changes in the direction 
or momentum affect the event classification. For exam- 
ple, TT^'s produced at the weak interaction vertex can be 
absorbed via intra-nuclear interactions within the tar- 
get nucleus, therefore escaping direct detection. Like- 
wise, intra-nuclear interactions can result in tt*^ produc- 
tion within the target nucleus, even in the absence of 
TT^'s at the weak interaction vertex. Therefore, the inter- 
actions of pions, kaons, etas and nucleons are also taken 
into account. The meson and nucleon interactions are 
treated using a cascade model, and each of the particles 
is traced in the nucleus until escaping from it. 

In our simulation, the following intra-nuclear pion in- 
teractions are considered: inelastic scattering, charge ex- 
change and absorption. The actual procedure to simulate 
these interactions is the following: first, the generated po- 
sition of the pion in the nucleus is set according to the 
Woods-Saxon nucleon density distribution [sl]. Then, 
the interaction mode is determined by using the calcu- 
lated mean free path of each interaction. To calculate 
these mean free paths, we adopt the model described by 
Salcedo et al. [531 ■ The calculated mean free paths de- 
pend not only on the momentum of the pion, but also 
on the position of the pion in the nucleus. If inelastic 
scattering or charge exchange occurs, the direction and 
momentum of the pion is determined by using the results 
of a phase shift analysis obtained from tt — N scattering 
experiments [H^. When calculating the pion scattering 
amplitude, the Pauli blocking effect is also taken into ac- 
count by requiring the nucleon momentum after the in- 
teraction to be larger than the Fermi surface momentum 
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at the interaction point. This pion interaction simulation 
is tested by comparison with data, including tt^^C scat- 
tering and pion photo-production + X) 
data. 

Re-interactions of the nucleons (protons, neutrons) 
produced in the neutrino interaction are also important. 
Each nucleon-nucleon interaction modifies the nucleon 
momentum and direction, possibly causing the number of 
visible nucleons to be mis-predicted if not properly mod- 
eled . Elastic scattering, single and two-pion produc- 
tion, are considered. 

Our simulation predicts that in 26% of SciBar CC in- 
teractions with TT*^ production at the primary neutrino- 
nucleon electroweak vertex, the neutral pion(s) does not 
escape the target nucleus. On the other hand, 15% of 
the events with tt^ emerging from the target nucleus are 
produced in nuclear interactions. 

B. Detector Response 

The GEANT3 ^ package is used to simulate the de- 
tector geometry and the interactions and tracking of par- 
ticles. The CALOR program library [6l| is used to sim- 
ulate the interactions of pions with the detector material 
for pions with momentum greater than 0.5 GeV/c. For 
lower energy pions, a custom library [4^ is used. 

The energy loss of a particle in each single SciBar strip 
and each individual EC sensitive fiber is simulated. The 
energy deposition is converted in the detector response 
taking into account the Birk's saturation of the scintil- 
lator, the light attenuation along the fibers, the Poisson 
fluctuation of the number of photo electrons, the PMT 
resolution, and the electronic noise. The crosstalk in 
nearby SciBar channels is also taken into account. 

In SciBar the timing of each hit is simulated from the 
true time of the corresponding energy deposition, cor- 
rected by the travel time of the light in the WLS fiber 
and smeared by the timing resolution. 

The MRD simulation includes both ionisation and mul- 
tiple scattering in the drift chambers. 

The input parameters of the detector simulation are 
derived from laboratory measurements and calibration 
data. The features of the simulation have been systemat- 
ically compared and tuned with cosmic ray and neutrino 
data. 



IV. EXPERIMENTAL SIGNATURE 

A. Definition of signal and background 

In this analysis the process we want to measure is de- 
fined inclusively with respect to a tt" produced in the 
target nucleus of the CC neutrino interaction. We call 
signal an event with one (or more) tt^ that comes from 
the neutrino interaction vertex or from a re-interaction 
inside the target nucleus. An event with an 77 decaying 



into one or more tt^ 's or into a 7 pair at the target nucleus 
is also considered a signal event. Events are considered 
background if the final state n° is only produced due to 
secondary interactions occurring outside of the target nu- 
cleus such as, for instance, charge exchange of a charged 
pion, or tt*' production in inelastic hadronic interactions. 
Another background category consists of the events se- 
lected accidentally where no tt*^ was produced. 

According to this definition, the CCtt" fraction pre- 
dicted by the neutrino Monte Carlo, integrated over the 
K2K energy spectrum, is 13.9% of the total number of 
neutrino CC interactions. The composition of the signal 
is the following: 

• 6.5% is resonant production: 5% with a tt*^ pro- 
duced in the resonance decay and 1.5% with tt" 
produced in a nuclear re-interaction in the target 
nucleus; 

• 6.6% is non resonant production, mainly DIS: 6.0% 
with one or more tt^'s produced at the neutrino in- 
teraction vertex and 0.6% in nuclear re-interactions 
in the target nucleus; 

• 0.8% comes from nuclear re-interactions, mostly 
CCQE, where a tt" is produced in the re- 
interactions in the target nucleus. 

The fraction of signal events with more than one tt" is 
43%. 

The final state topology of the CC inclusive events is 
characterized by one muon and at least two electromag- 
netic showers, plus possibly other particles coming from 
the neutrino interaction vertex. If the photon converts 
in SciBar, the hit patterns of the low energy electromag- 
netic showers are reconstructed by the SciBar tracking 
algorithm and the direction of the photon is given by 
the corresponding track. The SciBar conversion length 
is about 40-50 cm (SciBar in fact corresponds to 4 Aq). 
If the photon converts in the EC, the energy is recon- 
structed by the EC cluster algorithm and the position 
of the photon conversion is the energy-weighted average 
hit position in the cluster. Therefore the experimental 
signature is given by one track originating in SciBar and 
reaching MRD, at least two photons reconstructed ei- 
ther as SciBar tracks disconnected and pointing to the 
neutrino interaction vertex or as clusters in the EC. To 
isolate a sample of events that satisfy the topology de- 
scribed above, we first select a clean sample of CC events, 
characterized by a SciBar track matched with an MRD 
track. In Sec. IIVBI we describe the selection criteria to 
isolate the CC inclusive sample which is used for nor- 
malization. In the same section, we further classify CC 
events into sub-samples of varying CCQE purities, which 
are used to quote the CCtt" production cross section rel- 
ative to either the CCQE or the inelastic cross section. 
Out of the CC inclusive sample, we require further cuts 
to select photons and to separate the CCtt^ sample from 
other topologies. This is described in Sees. IIV CI and 
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IIVDI Section IIV El describes the reconstruction of the 
incoming neutrino energy. 



B. CC event selection 

The selection of a CC interaction requires a muon can- 
didate in the event. A muon candidate is a reconstructed 
3D track in the SciBar fiducial volume (FV) matching a 
reconstructed track in MRD. The FV is applied requiring 
the upstream edge of the track to be within ±135 cm in 
X and y, and -75 < z(cm) < 70 with respect to the center 
of SciBar. This corresponds to a 10.9 m'^ fiducial volume 
and 11.6 tons of fiducial mass. The track is also required 
to be in time with the neutrino beam, i.e. within ±50 ns 
with respect to the closest neutrino bunch. The extrapo- 
lation of the SciBar track is required to be matched with 
a track in MRD. The matching is with a MRD3D track 
or with a MRDIL hit, defined in Sec. Ill B 31 

The neutrino interaction vertex is reconstructed as the 
upstream edge of the muon candidate track in SciBar. 
The resolution in x and y is symmetric with 0.9 cm RMS. 
The resolution in z has a 1.6 cm RMS and a small satel- 
lite peak one SciBar layer (2.6 cm) upstream of the true 
neutrino vertex, due to crosstalk between MAPMT chan- 
nels. 

We select 11,606 events in the data and 432,856 in the 
full MC sample (before normalization) , with an estimated 
selection efficiency of 49.5% and a CC purity of 97.5%. 
The main background comes from the neutral current 
(NC) multi pion or single pion events in which a pion 
gives a signal in the MRD detector. The background 
induced from neutrons coming from the beam target is 
found to be negligible. 

In this analysis we consider four CC sub-samples, 
shown in Table IIIII which are characterized by differ- 
ent fractions of non-QE (nQE) and QE interactions. The 
first sample consists of events with a single reconstructed 
track and has 72.4% efficiency and 66% purity for QE 
events. For the events with two tracks, the direction of 
the second track is compared with the expected direc- 
tion of the proton in the assumption of a CC QE inter- 
action. If this angle A9p is smaller than 20°, the events 
are classified as "two tracks quasi-elastic". Events with 
A9p > 20° are further divided in two categories, depend- 
ing on whether the dE / dx of the second track is consis- 
tent with a pion or with a proton. The Monte Carlo is 
normalised to data using the first two samples in Table Hill 
which have the largest quasi-elastic contribution. The 
same normalisation is used in all plots before the fit. In 
all plots signal and different background components are 
stacked. In order to extract the result in Sec. |Vl we use 
the four samples described in Table IIIII and we leave the 
data to MC normalization free in the fit to properly ac- 
count for the correlation between the normalization and 
the other sources of systematic error. 



CC sub-sample 


<^QE% 


riQE% 


Data 


MC„ 


1 track 


72.4 


66 


6125 


6080 


2 tracks QE 


16.9 


76 


1262 


1307 


2 tracks nQE n 


2.3 


12 


1048 


960 


2 tracks nQE p 


7.3 


27 


1453 


1220 



TABLE III. Efficiency (eqb) relative to all QE events selected 
in the CC sample and purity (t?qb) for QE events, and number 
of events selected in MC and in data for the different QE and 
nQE samples 



C. Photon selection 

All CC selected events are subject to further selection 
criteria to tag photons. A photon candidate can be ei- 
ther a SciBar track or an EC energy cluster. In order 
to be considered as a photon candidate a SciBar track 
should satisfy the following requirements. First, the tim- 
ing of the track has to be within 10 ns with respect to 
the muon track; second, the track is required not to be 
matched with a MRD3D track. Third, the photon con- 
version point, defined as the track edge closest to the 
neutrino interaction vertex, is required to be within ±145 
cm in X and y and ±80 cm in z with respect to the SciBar 
center. Fourth, in both projections the distance between 
the photon conversion point and the neutrino vertex has 
to be larger than 20 cm; and fifth, the track extrapolation 
to the Z position of the neutrino vertex should be within 
25 cm from the neutrino vertex. The disconnection from 
the vertex of SciBar photon candidates, defined as the 
3D distance between the reconstructed neutrino vertex 
and photon conversion point, is shown in Fig. [5l The 
shape of the disconnection is consistent with MC. A fit 
with an exponential function between 50 cm and 130 cm 
gives a value of A^j^g^ = (51.8 ± 5.4) cm for data, in good 
agreement with A^jg^, = (52.9 ± 0.6) cm for MC. The re- 
sult is consistent with the electromagnetic origin of the 
selected photon sample. In fact, according to MC, 82% 
of the background events also contain a genuine photon. 

A fraction of photons converted in SciBar will have 
some energy leakage in the EC. Moreover, all the pho- 
tons not converted in SciBar and pointing to the EC will 
convert in the upstream (vertical) EC plane. We consider 
only clusters with energy larger than 50 McV for the ver- 
tical plane and 25 MeV for the horizontal plane. The en- 
ergy of the SciBar photon candidates and the associated 
EC vertical and horizontal clusters are added together in 
order to reconstruct the photon energy. EC clusters that 
do not match any of the reconstructed SciBar tracks are 
considered isolated. Isolated vertical clusters are paired 
to isolated horizontal clusters according to their energy 
and they are considered as additional photon candidates. 

Overall, 479 events with at least two photons are recon- 
structed in data and 380 in MC, with an overall efficiency 
of 7.6% and a purity of 59.2% (MC has been normalized 
to the data using the normalization factor described in 
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Ldisc (cm) 

FIG. 5. Disconnection: distance L^^jg^ between the recon- 
structed photon candidates conversion point and the neutrino 
interaction vertex. 



constructed in EC, and in 6% two or more photons are 
converted and reconstructed in EC. These relative frac- 
tions are weh reproduced by the MC. It is worth noting 
that this strongly supports the hypothesis that the excess 
of photon candidates is due to physics and not to detec- 
tion bias, since SciBar and EC are completely indepen- 
dent detectors, with different reconstruction efficiencies 
and systematics. 

As a crosscheck, we eye-scanned 100 data and 100 MC 
events. Despite the limited statistics and the subjectivity 
in the eye-scan classification criteria, the result is that the 
main features of the selected sample are well reproduced 
by the MC simulation. In particular, the background 
coming from the secondary interactions in data and MC 
agrees within the statistical uncertainties of this eye-scan 
cross-check. 



D. IT reconstruction 

The tt" mass is reconstructed from the energy and the 
direction of its two photon decay products: 



Sec. IIVB[) . Figure |6] shows the multiplicity of photon 
candidates per event. 
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FIG. 6. Photon multiplicity of CCtt" candidate events 



The excess of data with respect to the MC is 26 ± 6% 
(statistical error only). In 74% of the candidate events, 
all photons are converted and reconstructed in SciBar 
with possibly an energy leakage in EC. In 20% of the 
candidate events, one of the photons is converted and re- 



where E^-^ and E^^ are the reconstructed energies of the 
two photons, and ^71,72 is the opening angle between 
them. If the photon converts in SciBar, the direction is 
reconstructed using the SciBar 3D reconstructed track. 
If the photon converts in EC, we take as photon direc- 
tion the direction of the line connecting the reconstructed 
neutrino vertex and the center of the EC cluster. For the 
highest (lowest) energy photon reconstructed in SciBar, 
the energy and angular resolutions (FWHM) are 50 MeV 
(65 MeV) and 0.15 rad (0.18 rad), respectively. 

As shown in Fig. [HI in 21.4% of the selected events 
there are more than two photon candidates and therefore 
more than one tt" candidate. For these events the photon 
pair corresponding to the best 7r° candidate is selected as 
the combination which has the reconstructed tt" vertex 
closest to the neutrino interaction vertex. If there are one 
or more EC photon candidates (5.7% of the total sample) 
the best combination is selected as the photon pair with 
the reconstructed invariant mass closest to the tt" mass. 

Figure [7] shows the reconstructed tt" invariant mass for 
data and different MC contributions to signal and back- 
ground. The signal contribution (according to the defini- 
tion given in Sec. IIV Ap is divided into 7r° from resonant 
and non-resonant production and CCQE, and the back- 
ground is divided into resonant, non resonant production 
and CCQE plus neutral current (NC). It should be noted 
that most of the background contains a 7r° in the final 
state, so the shape of the invariant mass distribution for 
signal and background is similar. 
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tiq mass (MeV) 

FIG. 7. Reconstructed tt" mass before fit. 



E. Neutrino energy reconstruction 

The neutrino energy in a CC interaction can be recon- 
structed from the measured muon energy and angle using 
the following formula, provided the invariant mass W of 
the hadronic final state is known: 



j^rec ^ - O + '^E^{M,, -V)- (Af„ - Vf 

2x{-E^ + {Mn-V)+p^cos{9^)) ^' 

where V is the nuclear potential for carbon which is set 
to zero, and Pfj,, E^, and 0^ are the muon momentum, en- 
ergy and angle. For the QE final state, we have = 
and the formula used for neutrino energy reconstruction 
in the oscillation analyses. In the present analysis, 98% 
of the selected sample is non-QE, mostly resonant sin- 
gle pion production and DIS, and it is characterized by 
a broad W spectrum. We found W = 1.483 GeV the 
optimal value to reconstruct the neutrino energy in the 
MC sample of selected events. We use this value of W 
to reconstruct the neutrino energy in data and Monte 
Carlo. The uncertainties on the values assumed for W 
and for V will be considered as a source of systematic 
errors and evaluated in Sec. lVIl The resolution turns out 
to be 22%/^/E{GeV) for the selected sample. The as- 
sumption of an average W value is the largest effect in 
the reconstructed neutrino energy resolution. Using the 
true W value in Eq. [2]the resolution is 15%/ ^/E{GeV). 

The reconstructed neutrino energy is shown in Fig. [8] 
for data and different MC signal and background compo- 
nents. The threshold at about 1 GeV is due to the fixed 
value assumed for W in Eq. [51 




FIG. 8. Reconstructed neutrino energy before fit. 



V. LIKELIHOOD FIT 

From our sample of selected events we measured the 
ratio of the inclusive CCtt'^ cross section to the CCQE 
cross section. The uncertainties in the absolute neutrino 
flux cancel out in the ratio with respect to an indepen- 
dent and relatively well known process like the CCQE 
interaction. 

We performed a maximum likelihood fit of the recon- 
structed neutrino energy distribution of the CCtt*^ sample 
shown in Fig. [5] At the same time we fit the number of 
events in the different CC sub-samples described in Sec. 
IIVBI The likelihood function is defined as 

L = L^a ■ Lcc = WP{nk,p.k) • ]^P(ns,/^s) (3) 

k s 

where P{n, ji) is the Poisson probability for n observed 
events with expectation value /x. The maximum likeli- 
hood fit is calculated by minimizing the log-likelihood 
function F ~ —2Log{L) which follows a distribution. 

The index s labels the 4 CC sub-samples in Table Hill 
and the index k labels 50 bins spanning the range 0-5 
GeV of the reconstructed neutrino energy. 

The expected events fXs in each CC sub-sample (1- 
track, 2-track QE, 2-track nonQE pion and 2-track 
nonQE proton) are defined as: 

,, ft qQE , n ares , p nnres , 

+ Scc,J (4) 

where the number of Monte Carlo events contributing 
to signal (Sec) smd background (Bcc) a-re divided into 
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quasi elastic (QE), resonances production (res), non res- 
onant production (nres) and neutral current processes 
{NC). The non resonant production includes all CC non- 
QE processes different from resonant production, mainly 
deep inelastic scattering. 

The parameters Rres and Rnres are free in the fit in or- 
der to independently re-weight the corresponding Monte 
Carlo contributions relative to the quasi elastic process. 

All the MC distributions are normalized as described 
at the end of section IIVBI An additional overall normal- 
ization parameter /„ is left free in the fit. 

The number of expected events /i^ in Eq.[3]is given by: 



Fit Variable | Fit Result 
Energy independent fit 



RresSriEi)+RnresSr'%Ei) 



jyQE . 7j Tjres , p r>nres , rj 



NC 



(5) 



Sk and Bk are the Monte Carlo events, respectively sig- 
nal and background, contributing to the final CCtt^ sam- 
ple in each bin k of reconstructed neutrino energy. The 
Monte Carlo signal events are further divided in 4 bins 
El according to their true neutrino energy: 0-1.5 GeV, 
1.5-2.0 GeV, 2.0-2.5 GeV and greater than 2.5 GeV. The 
factors Fi are defined as: 



F, = 



Ek 


S^''{Et) + Sr{Et) + ST'^{El) 




Ek 




{Ei) + RresS^"^^ {Ei) -f RnresSJ^^*^^ 


{Ei)_ 



in order to keep the normalization of the signal events 
independent from Rres and Rnres- 

The fitting parameters are Rcc-K°{Ei) {j = 1,4), Rres, 
Rnres and /„. The best fit of Rcc-k°{EI) gives the dou- 
ble ratio data over Monte Carlo between the number of 
inclusive CCtt" events and the number of CCQE events, 
as a function of the true neutrino energy: 



R 



CCtt' 



'{Ei)^ 



N^SAEl)/N^§QE 



(6) 



The scaling of the inclusive CCtt^ contribution in the 
fit is energy dependent while the energy dependence of 
the CCQE is fixed to the Monte Carlo prediction since 
it has been accurately measured by previous experiments 
(.33 35, 33 )• The corresponding uncertainty is consid- 
ered a source of systematic error. We also performed 
an energy independent fit of the CCtt" to CCQE ratio, 
following the same approach as Eq. [S] but with a single 
fit parameter Rcctt" rescaling the CCtt'^ contribution re- 
gardless of the true neutrino energy. 





1.436±0.109 


Rres 


1.152±0.101 


Rnres 


1.373±0.241 


fn 


0.968±0.025 


Energy dependent fit 


RccTT"{Et) 


1.005±0.027 


RccTTo{Et) 


1.180±0.127 


RcCtto{Ev) 


1.307±0.198 


RccTTo{Et) 


1.418±0.129 


Rres 


1.105±0.098 


Rnres 


1.479±0.233 


U 


0.980±0.021 



TABLE IV. Energy dependent and independent fit results for 
Rcc.o{Ei). 



Table |TV] shows the best fit values of Rcc-k^ for 
the energy-independent fit and the the four parameters 
RccTT'>{Ei) for the energy-dependent fit. 

The xVd.o./. before the fit is 7135/44 = 162.1. The 
xVd.o./. for the best fit is 40.2/37 = 1.095 for the 
energy-dependent fit and 43.8/40 = 1.089 for the energy 
independent fit. 

The errors quoted for Rcctt" are purely statistical. 
The error induced on Rcctt'> by the absolute normal- 
ization /„ and by Rres and Rnres is evaluated in the fit 
in order to take in to account correlations but it is con- 
sidered a systematic error and reported in the first row 
of Table IVII together with the other sources of systematic 
in Sec. EB 

Figures |9] and \W\ show the reconstructed 7r° momen- 
tum and angle with respect to the beam direction in the 
laboratory frame, with the inclusive CCtt*^ production in 
the Monte Carlo rescaled to the best fit value for both 
signal and backgrounds. 

The fit results in Table IIVI show an excess of CCtt^ 
production with respect to our reference MC model. The 
energy dependent fit shows that the excess increases with 
the neutrino energy. The data to MC ratio for the non 
resonant processes Rnres is larger than one while the ra- 
tio for the resonant contribution Rres is consistent with 
one within the statistical uncertainty only. The reso- 
nant production with respect to the CCQE cross section 
was measured by the K2K collaboration in the CCl7r+ 
channel [13 and found to be: 0.7341°;^^°, m very good 
agreement with the MC prediction (0.740 ± 0.002(stat)). 
According to our reference MC model, 50% of the non 
resonant events have one or more tt'^ in the final state and 
44% of the selected CCtt^ sample is produced in non res- 
onant processes. Constraining the resonance production 
to the experimental value and uncertainty given above, 
we can use our CCtt^ sample to measure the non resonant 
contribution. 

We define CCnres (the CC non-resonant cross section) 
as the difference between the total CC cross section and 
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FIG. 9. 71° momentum distribution obtained rescaling the 
Monte Carlo witli tlie energy independent fit result. 
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FIG. 10. n'^ cos{9) distribution obtained rescaling the Monte 
Carlo with the energy independent fit result. 



Rnres 


1.461±0.118 


Rres 


1.112±0.098 


fn 


1.030±0.012 


Energy dependent fit 


Rnres ) 


1.064±0.266 


Rnres (-E'jy ) 


0.872±0.252 


Rnres (-£'1/ ) 


1.356±0.170 


Rnres (-S'jy ) 


1.567±0.164 


Rres 


1.109±0.102 


fn 


1.026±0.012 



TABLE V. Energy dependent and independent fit results for 

Rnres (^E^ ) . 



tion /„ were left free in the fit. The parameter Rres was 
also free and the experimental constraint was incorpo- 
rated in the fit by adding the term ^'^''^ji)^^'' 
Log- likelihood function derived from Eq.JS] 

Table |V] shows the best fit values for the energy inde- 
pendent fit of Rnres and for the energy dependent fit of 

Rnres {-^i/ ) ■ 

The xVd.o./. for the best fit is 41.8/39 = 1.07 for the 
energy dependent fit and 75.5/42 = 1.80 for the energy 
independent fit. The value of the x^/d.o./. for the en- 
ergy independent fit shows that a three parameters fit of 
signal and background, not taking into account the en- 
ergy dependence of the non-resonant contributions, gives 
a poor description of our data. 

This result is obtained assuming that tt^ and 7r+ pro- 
duction from resonances are constrained by the same pa- 
rameter within the Rein-Sehgal model. As a consistency 
check we repeated the fit without constraint to the CC7r+ 
measurement and the results are consistent within 3% 
with the results in Table [V] This difference is accounted 
by the systematic due to the uncertainty in the non-QE 
composition evaluated in Sec. I VII 

The errors quoted for Rnres are purely statistical. The 
error induced on Rnres by the normalization /„ and by 
Rres is evaluated in the fit in order to take in to account 
correlations but it is considered a systematic error and 
reported in the first row of Table IVTl The full systematic 
error is evaluated in Sec. IVII 



the sum of quasi elastic and resonance production. Using 
the CCtt" sample we perform an energy dependent and 
an energy independent fit of the CCnres to CCQE ratio, 
following the same approach as Eq.[5]but with Rcctt° and 
the normalization factors Fj fixed to one. The fit param- 
eter Rnres rcscaling the CC non resonant contribution 
in the energy independent fit (the parameters Rnres{El) 
in the energy dependent fit) and the overall normaliza- 



VI. SYSTEMATIC ERROR STUDY 

In this section we discuss the sources of systematic 
error. The contributors to the systematic error on the 
energy independent results CCtt^ and CCnres are sum- 
marized in Table rvTl The systematic errors for the energy 
dependent results in Tables IVlIl and IVIIII are calculated 
following the same approach. 
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Source 


crcCQE 


*^ CCnres tq^i 

[/o\ 

OCCQE 


Normalization and fit 


-1.8 -f 1.8 


-3.5 +3.5 


non-QE CC cross sections 


-3.4 H-3.3 


-3.1 +3.7 


Bodek&Yang corr. 


-4.3 -1-3.5 


-8.3 +7.8 


CCQE Ma 


-1.3 +2.4 


-0.8 +1.4 


NC/CC ratio 


-0.5 +0.5 


-0.8 +0.8 


u flux 


-0.1 +0.1 


-0.4 +0.4 


El, reco. parameters 


+0.2 -0.2 


+0.3 -0.3 


Int. model/flux 


-5.9 +5.7 


-9.6 +9.5 


n absorbtion 


-2.0 +2.1 


-1.8 +2.0 


n inelastic 


-3.0 +1.8 


-2.2 +1.5 


Proton rescattering 


-1.9 +0.3 


-2.8 +2.6 


Pion inter, length 


-1.5 +1.5 


-2.9 +2.3 


Nucl. model 


-4.1 +2.8 


-4.9 +4.3 


PMT resolution 


-0.5 +0.1 


-0.6 +0.7 


Scintillator quench. 


-0.1 +0.5 


-0.4 +0.5 


Cross-talk 


+1.2 +2.6 


+1.4 +2.3 


PMT threshold 


-1.7 +2.0 


-1.6 +2.2 


Detector effects 


-1.8 +3.3 


-2.2 +3.5 


Fiducial Volume 


-2.6 +2.5 


-3.3 +3.2 


V CI Lex OihC 






Vertex pointing 


-1.0 +1.8 


-1.4 +1.9 


EC cluster energy 


-0.4 +1.2 


-0.4 +1.1 


Selection cuts 


-3.4 +4.1 


-4.8 +4.7 


Total 


-8.1 +8.2 


-12.0 +12.0 



TABLE VI. Systematic errors for the CCtt" and the CC non 
resonant cross sections, relative to the CCQE cross section. 



A. Interaction model and neutrino flux 

The error coming from the absolute normalization for 
Rcctt" is estimated repeating the corresponding fit while 
fixing all the other parameters including the absolute nor- 
malization /„ at their best fit values. The resulting error 
for RccttO is the one reported in the previous section as 
pure statistical error while its quadratic difference with 
the full fit error is reported as normalization errors in 
Table IVIl The same procedure is applied to evaluate the 
errors on Rnres- 

The systematic error coming from the uncertainty in 
the composition of the non-QE CC cross section has been 
taken into account assigning a weight factor Wi to each 
non-QE CC channel and repeating the fit. A constraining 
function Fgyst was added to the Log-likelihood derived 
from Eq. O 



Fsyst — 



(Ty/2 



(7) 



where aw — 30%. The double ratio data to Monte Carlo 
of the NC to CC processes is also left free in the fit with a 
20% constraint, adding the corresponding term to Fgyst 
described in Eq. [T] The total systematic error listed in 
the second row of Table IVll takes into account the corre- 
lation between the different sources above. 

The QE axial mass is varied by ±10% (according to 
the uncertainties in the measurement reported in (62|) 



and the Bodek and Yang corrections to DIS events by 
±30%[43- The resulting systematic errors are added in 
quadrature in the total interaction model uncertainty. 

The uncertainty in the shape of the neutrino energy 
spectrum shown in Fig. [1] is considered by changing the 
fiux in each bin, taking into account their errors and the 
correlations between them (see Ref. 0] for details). 

The uncertainty on the values assumed for the neutrino 
energy reconstruction parameters in Eq. [5] is evaluated 
by changing the value of VF=1.483 GeV by ±15%, cor- 
responding to assume that all non-qe selected events are 
from A(1232) rather than from an average VF=1.483GeV 
The nuclear potential V is varied from to 27 McV. 



B. Nuclear model 

Nuclear effects alter the composition and kinematics of 
the particles produced in neutrino interactions in nuclei. 

Pion absorption and inelastic scattering processes, in 
particular pion charge exchange, modify the tt" yield. To 
account for the uncertainty in the Monte Carlo modeling 
of these effects, the pion absorption and pion inelastic 
scattering cross sections arc varied by ±30% [1^. The 
proton rescattering is changed by ±10% according to the 
uncertainties derived from cross section measurements 
[sl,!!!]- The systematic errors in Table IVll are calculated 
by repeating the analysis for each variation of the corre- 
sponding source. The uncertainty in the pion interaction 
length is considered by changing its value by ±20%. The 
overall uncertainty on the MC model is calculated con- 
sidering the uncertainty in the pion interaction length 
fully correlated to the pion inelastic cross section above. 



C. Detector effects 

The SciBar hit threshold, set at 2 photo-electrons, is 
changed by ±30% and the corresponding variation of the 
result is quoted as a systematic error. 

The model for the cross-talk in SciBar takes into 
account the second neighboring pixel and has a sin- 
gle free parameter n corresponding to the fraction of 
charge given by cross-talk in the adjacent pixel. The 
best fit obtained comparing data and Monte Carlo is 
n = (3.25 ± 0.01) X 10~^. The same model is used for 
the crosstalk simulation in Monte Carlo and for the cor- 
rection of the crosstalk effect both in Monte Carlo and 
data. To evaluate the systematic due to the crosstalk, 
we changed the crosstalk parameter n in the simulation 
in the range from 3.0% to 3.5%, corresponding to the 
uncertainty in the crosstalk modeling [l7| 

Smaller systematic detector effects are induced by the 
uncertainties in the single photo-electron PMT resolu- 
tion and the scintillator quenching (Birk's saturation) in 
SciBar. The SciBar PMTs resolution in the Monte Carlo 
is set at 40% [l^- This value was chosen by tuning the 
dE/dx per plane for muons in Monte Carlo to match the 
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response to cosmic ray data. The uncertainty is evaluated 
to be 10% and the corresponding systematic errors are 
hsted in Table IVII The scintillator quenching in SciBar 
was measured in a beam test and is well reproduced by 
Birk's eg nation [27}. The systematic error is evaluated 
by varying the Birk's parameter within its uncertainty. 
Other detector effects were found to give negligible con- 
tributions to the systematic error. 



Energy Range 
GeV 


Cross Section Ratio 

rr 1^ c 


> 0.0 


0.426 ± 0.032(stat.)± 0.035(syst.) 


0.0 - 1.5 


0.155 ± 0.039(stat.)±0.010(syst.) 


1.5 - 2.0 


0.577 ± 0.062(stat.)±0.037(syst.) 


2.0 - 2.5 


0.861 ± 0.130(stat.)±0.067(syst.) 


> 2.5 


1.627 ± 0.138(stat.)±0.103(syst.) 



D. Selection cuts 

In order to evaluate the systematic uncertainties due 
to the selection, in Table IVll we quote the dependence of 
the result to variations of the cuts. 

We change the fiducial volume by changing (simulta- 
neously in data and Monte Carlo) one at the time the 
fiducial volume cuts in the three coordinates according 
to the resolution for reconstructing the neutrino vertex: 
0.8 cm for both X and Y, and 1.6 cm for the Z. Then we 
add in quadrature the three corresponding variations of 
the result. 

The systematic uncertainty on the cut requiring the 
photon track to be disconnected from the vertex has been 
assessed by looking at the resolution on the neutrino ver- 
tex reconstruction and adding in quadrature the resolu- 
tion on the photon conversion point. We assumed the 
resolution on the photon conversion point to be equal to 
the resolution on the muon vertex. 

The cut on the photon track pointing to the vertex is 
applied to the distance by which the photon candidate 
track misses the vertex when extrapolated to the vertex 
plane. We take 0.14 rad for the 2D angle resolution of 
photon direction (0.12 rad and 0.14 rad, respectively, for 
the most energetic photon and the least one). The cut 
is applied on photons disconnected by more than 20 cm 
from the vertex. They have on average a 50 cm distance 
from the vertex. We set the variation as 50 • 0.14 = 7.0 
cm. In the cut region (25 cm), agreement between MC 
and data is quite satisfactory. 



VII. RESULTS 

Under the assumption that the detection efficiency are 
the same in data and Monte Carlo, the ratio between the 
inclusive CCtt^ cross section and the CCQE cross section 
can be calculated from Eq. [5] multiplying the best fit val- 
ues of Rcct^o given in Table ITVlbv the MC prediction for 
the cross section ratio in each neutrino bin E^: 



<7CCQE 



<^CCQE 



{ED 



(8) 



MC 



Table IVlIl shows the CCtt" cross section ratio to CCQE, 
integrated over all energies and as a function of the four 
neutrino energy bins. Figure [TT] shows the result as a 
function of neutrino energy. The vertical bars are the sta- 
tistical errors, the height of the filled areas corresponds to 



TABLE VII. Inclusive Cross Section Ratio 
tion of the neutrino energy. 



CCCQE 



as a func- 



Energy Range 
GeV 


Cross Section Ratio 

(ycCnres 
O-CCQE 


> 0.0 


0.419 ± 0.034(stat.)± 0.050(syst.) 


0.0 - 1.5 


0.010 ± 0.002(stat.)±0.002(syst.) 


1.5 - 2.0 


0.432 ± 0.125(stat.)±0.056(syst.) 


2.0 - 2.5 


1.304 ± 0.164(stat.)±0.117(syst.) 


> 2.5 


2.954 ± 0.309(stat.)±0.354(syst.) 



TABLE VIII. Cross Section Ratio '^'"^ '^^ as a function of 



the neutrino energy. 



lyccQE 



the statistical and systematic errors added in quadrature 
and the data points in each bin are set at the weighted av- 
erages of the true neutrino energy for the selected CCtt" 
events. Figure [Til also shows the CCtt'^ over CCQE ratio 
and the two largest contributions, single pions from reso- 
nances and pions produced in DIS, as they are predicted 
by our reference MC. According to our reference MC 
the average true neutrino energy for the selected CCQE 
events is I.l GeV. The average true neutrino energy for 
our selected tt" sample is 1.3 GeV and 2.5 GcV for the 
fraction of tt" produced in DIS events. 

The CC non resonant to CCQE ratio is obtained from 
the best fit values of Rnres in Table |V] similarly to Eq. (5) 



<^CCr. 



(^CCQE 



-{ED = R^resiEi) X 



(^CCQE 



-{ED 



(9) 



MC 



The results for are reported in Table I Villi 

^CCQE ^ 

integrated over all energies and as a function of the four 
neutrino energy bins. 



A. Comparison with other experiments 

Past experimental results exist for the exclusive Vyji — > 

M~ TT^ cross sections on deuterium (Barish , Radecky 
and Kitagaki [l3l)- There is also a published result 
for the exclusive cross section v + p ^"pir'^i:^ (D^'Y 
[31). In order to compare with our result on CgHg, cross 
sections on deuterium have been rescaled to the differ- 
ent number of protons and neutrons. The ratio between 
CCtt^ and CCQE cross sections is computed by divid- 
ing the experimental results quoted above by the CCQE 
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cross section measured by Barish 



Below 1.5 GeV 



neutrino energy our result can be directly compared with 
the published single pion cross sections, since this is the 
main contribution to the inclusive cross section. The 
three points shown as diamond-shaped symbols at 1.07 
GeV, 1.70 GeV and 3.0 GeV are obtained adding the two 
pion urpn~^TT^ from [l^ to the single pious ji~pTr^ taken 
from [ij]. 
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O 



O 
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O 

1.2 
^ 1 
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K2K MC single pion 
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1 2 3 

Neutrino energy CCQE(GeV) 



FIG. 11. CCtt" to CCQE cross section ratio as 
a function of neutrino energy. The result of this 
analysis (•) is compared with our standard MC ex- 
pectation and past experimental results. The previ- 
ous experimental data are: (Q]ANL82[l3] ,(□) ANL[33 , 
(A)BNL86[l2,(0)ANL83-fBNL86[i||[l2(see text for de- 
tails). 



The MiniBooNE collaboration has recently published 
charged current tt*' production [6^ at a slightly lower 
neutrino energy (between 0.5 GeV and 2 GeV) and their 
result is consistent with the one presented in this paper. 



VIII. 



CONCLUSION 



Out of a sample of 11,606 charged-current neutrino in- 
teractions in the K2K SciBar detector, we selected 479 
inclusive tt^ events with an efficiency of 7.8% and a purity 
of 66.5%. The sample corresponds to 2.02 x 10^^ protons 
on target recorded with the SciBar -l-EC-l-MRD detectors 
at 1.3 GeV average beam neutrino energy. We measured 
the cross section for charged-current tt*^ production rela- 



tive to the charged-current quasi-elastic cross section to 
avoid the large uncertainties in the absolute neutrino flux 
determination. The result integrated over the neutrino 
energy spectrum is 



f^CCQE 



0.426 ± 0m2{stat) ± 0.035(sysi.) 



higher than the prediction of our reference Monte Carlo. 
The energy dependent CCtt'^ to CGQE cross section ra- 
tio is presented in Tabic IVIII and shown in Fig. [TT] The 
results of the best fit for the composition of our CCtt" 
sample show that the data excess comes from non reso- 
nant processes, mainly tt*^ production in DIS, rather than 
from tt'^ in resonance production. Using the measured 
CG single charged pion cross section [l3| as a constraint 
for the resonant production, we measured the ratio be- 
tween the CC non resonant and the CCQE cross section, 
integrated over the neutrino energy spectrum: 



C^CCQE 



0.419 ± 0.034(siai.) ± 0.050{syst.) 



For CC non resonant processes we define any charged 
current process except quasi elastic interaction and reso- 
nances production. The energy dependent cross section 
ratio is presented in Table I Villi The results presented 
here are the firsts for neutrinos of few-GeV energy on a 
GgHg target material and improve the precision of pre- 
vious results on different targets and therefore are a sig- 
nificant contribution to the knowledge of neutrino inter- 
action processes relevant for several present and future 
oscillation experiments. 
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